Shot noise in electron transport through a double quantum dot: A master equation 

approach 
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We study shot noise in tunneling current through a double quantum dot connected to two electric 
leads. We derive two master equations in the occupation-state basis and the eigenstate basis to 
describe the electron dynamics. The approach based on the occupation-state basis, despite widely 
used in many previous studies, is valid only when the interdot coupling strength is much smaller than 
the energy difference between the two dots. In contrast, the calculations using the eigenstate basis 
are valid for an arbitrary interdot coupling. We show that the master equation in the occupation- 
state basis includes only the low-order terms with respect to the interdot coupling compared with 
the more accurate master equation in the eigenstate basis. Using realistic model parameters, we 
demonstrate that the predicted currents and shot-noise properties from the two approaches are 
significantly different when the interdot coupling is not small. Furthermore, properties of the shot 
noise predicted using the eigenstate basis successfully reproduce qualitative features found in a 
recent experiment. 
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I. INTRODUCTION 

Precise control of coherent coupling between quantum 
states is of great importance in quantum information pro- 
cessing. Recent studies show that artificial two-level sys- 
tems designed using mcsoscopic circuits can be controlled 
in nanosecond time scales and can also exhibit coher- 
ent oscillations between two quantum states (see, e.g., 
Rcfs. |l|i2iji^). A double quantum dot (DQD) provides a 
useful system to explore coherent effects because inter- 
dot hopping intrinsically couples states in two different 
dots and is tunable via the gate voltage4i^ A commonly 
used observable for studying the effects of coherent cou- 
pling is the current through the DQD. Recently, shot- 
noise measurement has recently been demonstrated as 
another useful tool to study the coherent effectsi^ii More- 
over, the shot-noise properties have been predicted to be 
an indicator of the degree of entanglement between elec- 
tron states^i^ and they are also related to the radiative 
decay properties of the one-demensional quantum ring 
excitoniiS 

Shot noise, i.e., current fluctuations due to the dis- 
crete and stochastic nature of electron transport, de- 
scribes the correlation between electrons transported suc- 
cessively through mesoscopic systems, such as quan- 
tum dots (QDs) or molecular devices (for reviews, see 
Refs. [TTIIT2h . In classical transport, the noise is typi- 
cally Poissonian with a power density S = 2e(/), where 
e is the unit charge and (/) is the average current. 
However, either Coulomb interaction or the Pauli's ex- 
clusion principle can induce a negative correlation be- 
tween successive transport events. This reduces the noise 
power density so that S < 2e{I) corresponding to a sub- 
Poissonian noise.— In contrast, the interplay between 
Coulomb interaction and the Pauli's exclusion principle 
can also produce a positive correlation between the trans- 



port events, i.e., S > 2e(/). This corresponds to a super- 
Poissonian noise. The Fano factor F = S/2e{I) is usu- 
ally used to characterize the shot noise, where F ~ 1, 
F > 1, or F < 1, respectively, corresponds to the Poisso- 
nian, super-Poissonian, or sub-Poissonian noise. Many 
theoretical works show that super-Poissonian noise of 
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and positive cross cor- 



relation between different spin states^ in QDs can be 
induced via dynamical channel blockade. Moreover, a 
super-Poissonian noise in tunneling current caused by dy- 
namical channel blockade has been observed in a system 
consisting of two electrostatically coupled QDa^^ and also 
a single Q d23,24 -^^ recent experiments. 

In this work, we study the current and shot-noise in 
electrons tunneling through a DQD. We apply two dif- 
ferent approaches and compare the results. First, we fol- 
low many previous investigations (see, e.g., Refs.QIHIH^ 
and derive a master equation for the electron transport 
based on the occupation-state basis of the DQD. We 
show that to arrive at this master equation, one needs 
to assume that the interdot coupling strength is much 
smaller than the energy difference between the two dots. 
However, using realistic model parameters for the DQD, 
only Poissonian or sub-poissonian shot noise is predicted, 
while super-poissonian noise was also observed in a recent 
experimenti^ 

Alternatively, we derive a more generally applicable 
master equation in the eigenstate basis of the DQD, 
which does not require the assumption of a small interdot 
coupling and is hence valid for any arbitrary interdot cou- 
pling strength. The two master equations are formally 
different in general and are identical only in the limiting 
case when the interdot coupling is much smaller than the 
energy difference between the two dots. We show that for 
small interdot coupling, the properties of the shot noise 
predicted by the two master equation agree with each 



other as expected. However, for large interdot coupling, 
they are significantly different. More importantly, for 
typical model parameters, the shot noise deduced using 
the master equation in the eigenstate basis exhibits rich 
properties including Poissonian, sub-poissonian as well 
as super-poissonian statistics in good agreement with re- 
cent experimental observations (Ref. ly). Furthermore, 
qualitative features of the current and the shot-noise can 
easily be explained intuitively using the master equation 
in the eigenstate basis. 

The present paper is organized as follows. In Sec. II, 
we introduce the model for a DQD connected to two elec- 
tric leads. A phonon bath that affects the dynamics of 
the DQD is also considered. Two master equations for 
the electron dynamics in the DQD are derived in both 
the occupation-state basis and the eigenstate basis. In 
Sec. Ill and Sec. IV, we study, respectively, the proper- 
ties of the current through the DQD and the associated 
shot noise. Results based on the two master equations 
are compared. In Sec. V, we discuss the relation be- 
tween the two master equations. A brief conclusion is 
presented in Sec. VI. Finally, Appendixes A and B give 
detailed derivations of the master equations in both the 
occupation-state basis and the eigenstate basis. 



II. TIME EVOLUTION OF THE REDUCED 
DENSITY MATRIX OF A DOUBLE QUANTUM 
DOT 

The schematic diagram of a DQD connected to two 
electrodes by tunneling barriers is shown in Fig. [TJ The 
voltage across the DQD is biased so that the chemical 
potential of the left electrode /ii is higher than that of 
the right electrode ^n. Thus, an electron can tunnel 
from the left electrode to the right one via the DQD. We 
assume that the DQD is in the Coulomb regime (with 
both strong intra- and interdot Coulomb repulsions), so 
that at most a single electron is allowed in the DQD. In 
the occupation representation, the electron basis states 
are the vacuum state |0), the state with one electron in 
the left dot |1), and the state with one electron in the 
right dot 1 2). 

The Hamiltonian of the whole system reads (taking 
h = l) 





FIG. 1: (Color online) Schematic diagram of an electron 
transported through a DQD connected to two electric leads 
via tunneling barriers, (a) Considering DQD electron states 
in the occupation-state basis, the electron tunnels sequentially 
from the left lead to the right lead via first the left dot and 
then the right dot. (b) Considering the eigenstate basis, the 
electron is transport via either the ground-state channel or 
and the excited-state channel. The efi'ective tunneling rates 
from the left lead to the ground state and the excited state 
are FlQ^ and TlP^, while those from the ground state and 
the excited state to the right lead are r_R/3^ and Tro^. 



(a = l,r); 



ala2 



a[ai and ax 



a|a2 



are the Pauli matrices, with a[ (aj) being the electron 
creation operator in the left (right) dot of the DQD. In 
Eq. ([U, the first term saz/2, with e = 62 — £i denot- 
ing the energy difference between the two dots, gives the 
Hamiltonian of the two uncoupled quantum dots, while 
the second term fla^ characterizes the interdot hopping. 
The tunneling coupling between the DQD and the elec- 
trodes is described by 



ii^tot =i^lcads + HuQD +Ht + Hpi, + Hep. (1) Ht = ^{fllk a\ Cik + flrk 4 '^r Crk + H.C.) 



(4) 



The first two terms i^loads and -ffDQD arc respectively the 
Hamiltonians of the two electrodes and the DQD, and are 
given by 



where 



lk(rk) 



is the tunneling strength between the QD 



ak 

^^DQD ~ -Z^z+^(^x, 



(2) 
(3) 



where c^^^^ (Caka) is the creation (annihilation) opera- 
tor of an electron with momentum k in the electrode a 



and the left (right) electrode. The operators (Tj!) de- 
creases (increases) the number of electrons having tun- 
neled into the right lead (via the barrier between the 
DQD and the right lead).— These counting operators al- 
low one to keep track of the tunneling process during 
the evolution of the DQD. Below we focus on current 
and shot noise in electron tunneling through the right 
tunneling barrier, so only the related counting operators 
(Tr and Tj.) are introduced in the tunneling Hamiltonian 
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Also, we consider the effects of the phonon-bath envi- 
ronment on the evolution of the DQD. The Hamiltonian 
of this phonon bath is 

9 

with &J (bq) creating (annihilating) a phonon with fre- 
quency LOq . The electron-phonon interaction is given by 

H,p = a,Y,Xq{bl + bq), (6) 
q 

where Xq is the electron-phonon coupling strength. 

The evolution of the whole system is described by the 
von Neumann equation for the density matrix pj^ of the 
whole system: 

pjiit) = -i[Htou PR{t)]- (7) 

Here we are interested in the time evolution of the DQD 
and treat both the electric leads and the phonon bath as 
the total outside environment. We will hence derive the 
master equation of the reduced density matrix pd of the 
DQD: = Tr^l/Oi^}, where Tr^j - • • } denotes the trace 
over the degrees of freedom of both the electric leads 
and the phonon bath. In our calculations, we adopt the 
interaction picture based on the free Hamiltonian 

Ho = -ff leads + -f^ph + -f^DQD, (8) 

and the interaction Hamiltonian becomes 

Hint{t) - HT{t)+H,p{t), (9) 

where operators in the interaction and the Schrodinger 
pictures are related by A{t) — ^iHat j^^-iHat 
erator A. 

After tracing over the degrees of freedom of both the 
electrodes and the phonon bath, one obtains the master 
equation for the reduced density matrix of the DQD 
in the interaction picture as2^ 

= -*Tr£[ifi„t(t), Pd{0)pEm 

JQ 

(10) 

where p^; is the density matrix of the outside environ- 
ment. Because the trace of a single unpaired creation or 
annihilation operator over the lead or the phonon bath 
is zero, e.g., TTElcikPs} = 0, the first term in Eq. PU]) 
vanishes. Within the Born-Markov approximation, we 
have 

p'dit) = -Tl-E f dt'[H-Ut), [Hintit'), p'd{t)pEm ], 

Jo 

(11) 



Here the Born approximation amounts to the use of the 
second-order perturbation theory with respect to the in- 
teraction Hamiltonian iJint, while the Markov approx- 
imation assumes that the correlation times of the out- 
side environment (both the electric leads and the phonon 
bath) are much shorter than the typical quantum-state 
evolution time of the DQD. 

Since the left lead, the right lead and the phonon bath 
are completely independent of each other, the density 
matrix of the outside environment pE can be written as 
a tensor product of density matrices that describe the 
subsystems, i.e., ps = PLPRPph, where pL, Pr and pph 
are, respectively, the density matrices of the left lead, the 
right lead and the phonon bath. Therefore, the trace of 
the integrand in Eq. (fTT|) can be expressed as: 

Tr£[i7i„t(i), [i?i„t(i'), P^Wpb] ] 
= J2 Tia[i/T(0, [H^{t'),p'd{t)pa]] 

OL—l,r 

-f Trph [i?ep(i), [Hepit'), P^(i)Pph] ]. (12) 

Equation pip can thus be written as a sum of two cor- 
responding parts: 

piit) = CTp'At) + Cpi^piit). (13) 

Here, the dissipativc part due to the electric leads is given 

by 

t 

^TPd(0-- T^^'c. J dt'[HT{t)HTit')p'd{t) Pleads 

a=l,r 

-HTit)p'dit) Plcadsi?T(t') + H.C. ] , (14) 

where pioads = PlPr is the density matrix of the two 
electric leads. The dissipativc part caused by the phonon 
bath is 

t 

>CphPd(0 =-Trphy di'[ffep(i)i?cp(i')Pd(t)Pph 



~H,p{t)p'^{t) pphi?cp(t') + H.C. ] . (15) 

From Eqs. pS jl -lfTS ]) . one can derive the master equa- 
tion for the n-resolved reduced density matrix p^"^ (t) of 

the DQD, where plP\t) = {n\pd{t)\n) and n is the num- 
ber of electrons that have arrived at the right lead at time 
t. Below we derive two versions of the master equation in 
both the occupation-state basis and the eigenstate basis 
and then use them independently to study the current 
and shot-noise properties of the DQD. 

A. Master equation in the occupation-state basis 

The master equation of the DQD in the occupation- 
state basis was previously used to study the current 
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propertics^i2& and shot-noisc propcrtie o^'^^i"'^^'^" of elec- 
trons tunneling through the DQD. The occupation-state 
basis is defined by the states |0), |1), and |2), which corre- 
spond to the states of an empty DQD, one electron in the 
left dot, and one electron in the right dot, respectively. 
In the interaction picture defined by the free Hamilto- 
nian Hq in Eq. ([8|), the unperturbed evolution operator 



Uoit) 



is difficult to calculate in the occupation- 



state basis in the presence of intcrdot coupling, 
hence split into two parts: 



Ho — Hi + Hi', 



where 



Hi — iJlcads 

Hn = Vldx- 



Hph 



One 
(16) 

(17) 
(18) 



Following previously work o^^'^° on deriving master equa- 
tions in the occupation-state basis, we assume that the 
interdot couping Q is small and satisfies <C |e2 — ei|, 
we have Hq ~ Hi. The evolution operator can then be 
approximated as 



JHit 



(19) 



With this approximation, one easily obtains 

7?T(t)«5][f^a-alQfce-'("--"^)* 
fc 

-f ^rk T,. 4 c,fce-'('"'-'=-"^)* + H.c. ], (20) 

and 

H,p{t) « {a+e''' + a_e''-^*) ^ A,(6te'-''* + fo^e^'^"'*), 

(21) 

where 0Ji,2 = T£/2. Here aj^ ~ a\ai and (t_ = cr^ arc 
the raising and lowing operators in the occupation-state 
basis. 

We now consider the nonequilibrium case with a large 
bias voltage across the DQD, so that all energy lev- 
els of the DQD lie within the bias window, as shown 
in Figs, ma) and [Hb). Substituting Eqs. ^ and 
PT|) into Eqs. P^ - ([TC)) . taking the trace over the de- 
grees of freedom of both the two electrodes and the 
phonon bath, and converting the obtained equation to 
the Schrodinger picture, one arrives at the master equa- 
tion in the occupation-state basis for a weak interdot cou- 
pling ^ |e2 — Ell (see Appendix A): 

Pd{t) = -«[^DQD, Pd] + —'D[a\]pd + —'D[Tla2]pd 



(22) 



where ^l(r) = '^''^Pik{rk)^fk(rk) electron tunneling 

rate through the left (right) tunneling barrier. Here, the 



electron density of states pak at lead a {a = l,r) and 
the tunneling strength flak are assumed to be energy- 
independent. The notation 2? acting on any operator A 
is defined as 



V[A]p = 2ApA'^ - [A^Ap + pA^A]. 



(23) 



The dissipation rates induced by the electron-phonon in- 
teraction are 



71 

72 



27r{ J (e) [n (e) + 1] + J {-e) n (-e)}, 
2n{j i-e) [71 i-e) + I] + J (e) n (e)}, (24) 



where 



q 



(25) 



is the bath spectral density and n (w) = 
[exp (w/fc^T) — 1]~^ is the average phonon number 
at temperature T. Using Eq. p2|) and the relations!^ 



{n\'TlTrPd\n) 
{n\'TrTlpd\n) 
{n\TlpdTr\n) 
{n\TrPd^l\n) 



Pd ' 

(n) 
Pd ' 

(n-1) 
Pd > 

(n+1) 
Pd ' 



(26) 



one obtains the equation of motion for each density ma- 
trix element: 



Poo (0 ~^LPoa + ^RP22 

Mi^ {t) = T,p^^+^n{p^-^-pt;> 



p^^ {t)^-Tnpt^^7n[p^!^ 



(n) 
P2I 



(n) (n) 

71P22 -72Pii , 

(ti) (n) 
- 71P22 +72P11 , 



p';^'{t)^^ep^;$+^n[pTl'-p'^']- 



(ti) (n)\ ^R+ll+ 72 (n) 



P12 ■ 

(27) 



Then, the ith diagonal matrix element pu = ^„ p^"^ {i = 
0, 1, or 2) gives the occupation probability of the state 

The off-diagonal matrix element pi2{t) = P12 
describes the coherence between states |1) and |2), and 
P2i{t) = Pi2W- This master equation was used in many 
previous studies, e.g., Refs. I25I andl26l. 

The physical meaning of the master equation can be 

(n) 

understood as follows. Take the equation for p^i in 
Eq. (p7)) for example. The first term on the right-hand 
side describes the process of an electron tunneling from 
the left lead to the left dot with rate T^. The second 
term represents the coherent coupling between states 1 1) 
and 1 2) due to the interdot coupling. The third term de- 
scribes the phonon-induced relaxation process from state 
1 2) to |1) with rate 71. Finally, the fourth term describes 
the inverse process with rate 72. 



5 



B. Master equation in the eigenstate basis 

As explained above, the master equation in the 
occupation-state basis is only valid for a weak interdot 
coupling. To extend the results to an arbitrary interdot 
coupling f2, we now derive the master equation in the 
eigenstate basis of the DQD. The result is valid for any 
arbitrary interdot coupling strength. 

Diagonalizing the Hamiltonian of the DQD [Eq. ([3])], 
one has 



^DQD = ^(|e>(e|-|.g>(.g|) = %a(^\ 



2 



(28) 



where ilg — \/e^ + 4il'^ is the energy splitting of the two 
eigenstates of the DQD given by 



|e) = sin-|l)+cos-|2), 



l.9> = cos-|l)-sin-|2), 



(29) 



with tan6' = 2f7/e. The eigenstates and the occupation 
states are related by 



|1) = cos - \g) + sin - |e) , 

6 
|2) = -sin - |.g) +C0S- |e) 



(30) 



With these relations, the tunneling Hamiltonian [Eq. ([4])] 
and the electron-phonon interaction [Eq. ([6])] can be writ- 
ten, in the eigenstate basis, as 



^ik ( cos 2 + 2 ) 



■ sin 2 «g + cos -j^a\)Tr Crk + H.c. 



sin 9 ai"^ + cos a^^A \ {^1 + 



(31) 

In the interaction picture based on the free Hamiltonian 
Hf) given by Eq. they become 



HT{t) = E ( cos ^ 4 e»-«* + sin ^ fltg-e* ' 



Clk 



• Qrk [- sin -ate'"»* + cos - ale"^'* ' 



xT, 



*+H.c.}, (32) 
,{t) = [sin0 tx^^' +cos^ (a^^'e*"«* -Haf!^e-'"°*) 



bae~ 



(33) 



Now, one can evaluate Eqs. ([T3)) - (|T5l) using Eqs. ((321) 
and ([55)1 . After converting the result to the Schrodinger 
picture, one obtains the master equation in the eigenstate 
basis which holds for any arbitrary interdot coupling £7 
(sec Appendix B): 



Pd{t) = -i[i?DQD, Pd{t)] 



Ai, 



(34) 



where 



a = cos - 



2rjn 



^ = sin - 



2r!o 



(35) 



and 



Ai = 7ocos2 6'[n(f7o) + l], 
A2 = 70 cos^ 0n(rJo)- 



(36) 



with 7o = 27rJ(rio). Using Eq. ((M)) . the n-rcsolvcd equa- 
tion of motion for each density matrix element can be 
written as 



p^:^{t) = VlP^p^^ - Tro^p^^ - Aip(;') + A24";. 



(37) 

It follows from Eq. ([571) that the effective tunneling rate 
from the left lead to the ground (excited) state \g) (|e)) is 
Tlo^ (Fl/?^), while the effective tunneling rate from the 
ground (excited) state to the right lead is Tii0^ {Tb.cP') 
[see Fig. 1(b)]. We emphasize that these results derived 
in the eigenstate basis are valid for any arbitrary interdot 
coupling. 



CURRENT THROUGH THE DOUBLE 
QUANTUM DOT 



III. 



To compare the master equations in Eqs. ((27|) and ((37|) 
derived, respectively, in the occupation-state basis and 
the eigenstate basis, we first apply them to study the 
tunneling current through the DQD. In the next section, 
the associated shot noise will also be studied. The current 
I{t) through the DQD at time t is given by 



m 



dN{t) 
dt 



(38) 



where Wg.e = =Ff7o/2. Here cri^^' — agO^ and cr'j^^'' ~ (c'^'')^ 
are the lowering and raising operators in the eigenstate 
basis. 



where N{t) is the number of electrons that have tunneled 
into the right lead. Here, i is summed over all basis states 
of the basis used. 



6 




where 



0.00 



-200 200 

8j - 8^ (|ieV) 



FIG. 2: (Color online) Stationary current I" through the 
DQD as a function of the energy difference ei — £2 calculated 
using the occupation-state basis (7" = lo) and the eigenstate 
basis (I" ~ Ic) for a large interdot coupling ^ — 32 fieV. We 
have taken Tl ~ 100 ^eV, Fr = 2.5 fieV, 70 = 0.6 ^eV, and 
T = 2 K. 



Denoting results based on the occupation-state basis 
and the eigenstate basis by "o" and "e", values of the 
current Io{t) and Ic{t) calculated using Eqs. (|27|) and 
([57)1 are 



Io{t) = eTRp22{t), (39a) 
Ic(t)^e[ rR(3^pgg {t) + Tna^p^e (t) ] , (39b) 



At steady-state virith pii{t) = 0, calculated values and 
of the stationary current are 



eTLTi 



X {72 [46^ + (71 + 72 + rfi)2] 



A 

+ 4^(2(71 +72+ r^J}, (40a) 
eTLTn[P'Xi+f{X2+P'rR)] ^ ^^^^^ 



A= 4rj2r^ (2rL + r^) + (72 + r^) Tr (ie^ + r?^) 



+ (71 + 72) 



4rL(e2 + 2f72) + 4r22rfl 



+ ^lIii + l-if + ^r{12 + 3ri)(7i + 72) 



r|(272 + 3ri) 



S= A2 (Fl + a^Tn) + Ai (Fl + P^Tr) 



(41a) 



(41b) 

Figure [5] shows the calculated values and of the 
stationary current through the DQD. We choose a typical 
interdot coupling = 32 /ieV, which is experimentally 
accessible For both and at the resonant tun- 
neling point characterized by £1—62 = 0, the current 
reaches its maximum. Moreover, it can be seen that the 
current is asymmetric around the maximum point. This 
asymmetry was also observed in a recent experiment by 
Barthold et al.^ It is due to dissipations induced by the 
phonon bath, as we will now demonstrate. In the absence 
of electron-phonon coupling, we have 71 = 72 = and 
Ai = A2 = 0, so that Eqs. (|40ap and (|40bp reduces to 



4172 (2rL + Tr) + ' 
ien^TLTR 



(42a) 



(42b) 



^ 4172 (2ri + Tr) + 4e2FL ' 

where Eq. (|42ap agrees with the result from previous 
studieSfii^ in which the occupation-state basis was also 
used. From Eqs. (|42a[) and (j42bp . it is clear that the cur- 
rent as predicted using either the occupation-state basis 
or the eigenstate basis is symmetric about the current 
peak at £1 = £2- This reveals that this asymmetry of the 
current is due to the coupling of the DQD to the phonon 
bath. 

As shown in Sec. II, the master equation derived in the 
occupation-state basis is only valid in the limit of a weak 
interdot coupling, i.e., fl ^ |£2 — £i|, while that in the 
eigenstate basis is valid for any arbitrary interdot cou- 
pling. Figure [3] plots the values of the stationary current 
calculated in the two bases for a small interdot coupling 
(fi = 1 ^eV). As expected, when the interdot coupling 

is much smaller than the energy difference |£2 — £i| 
of the two dots, the stationary-state current calculated 
in the occupation-state basis agrees very well with that 
in the eigenstate basis, as is evident from Fig. [31 How- 
ever, when the interdot coupling is comparable to the en- 
ergy difference, the stationary currents in the occupation- 
state basis deviates drastically from the stationary cur- 
rent in the eigenstate basis (see Fig. [5]). This deviation 
can also be revealed in Fig. [3] in the narrow region with 
|£2 - £i| (= lp,eV) (see the inset of Fig. [S]). These 
clearly show the inaccuracy of the current and hence the 
master equation in the occupation-state basis at large Q. 
In this case, one must use the master equation derived in 
the eigenstate basis. 
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FIG. 3; (Color online) Stationary current /" as a function 
of £1 — £2 calculated using the occupation-state basis [P — 
Jo) and the eigenstate basis {I" = /|) for a small interdot 
coupling f2 = 1 /leV. Inset: The enlarged diagram of the 
stationary current P in the region with \e2 — £i \ comparable 

to n. 



IV. SHOT NOISE 

To calculate the shot noise in the tunncUng current 
through the DQD, it is particularly useful to define a 
generatin g fu nction for an electron counting variable s 
(see Refs.lllilandlll): 



G(t,5) = ^/y")(t). 



(43) 



and the variance reads 



dHrG{t, 1) 



+ (n) - (47) 



Applying the Laplace transform to the equation of mo- 
tion, Eq. ((44|). of the generating function, one has 



(48) 



Because of the incoherent long-time stability of the con- 
sidered system, the real parts of all the non-zero poles 
of G{z, s) are negative. Therefore, the long-time be- 
havior is determined by the pole zq closest to zero, i.e., 
G{t, s)^ g{s)e^"* . By the Taylor expansion of the pole 



rn>0 



(49) 



one obtains 



Clt, 



Mil 

ds 

d^g(l) /%(!) 
ds'^ \ ds 



+ (ci + 2c2)t. 



(50) 



In particular, the Fano factor of the shot noise is given 
by 



This generation function obeys the equation of motion 

G{t,s) ^ M{s)G{t,s), (44) 

where M{s) is a transition matrix that can be calculated 
using the master equation [Eq. or Statistics on 

the number of transported electrons n can be determined 
from the derivatives of the generating function: 



ffltrG(t, 1) 



In particular, the mean of n is 

dti-G{t, 1) 



in) 



ds 



(45) 



(46) 



F=l + 2^, 
ci 



(51) 



where F > 1 {F < 1) indicates super (sub)-Poissonian 
noise, compared to = 1 for classical Poissonian noise. 

We first consider results based on the occupation-state 
basis. To calculate the Fano factor, we can show, using 
Eqs. (H?]) and (gS]), that the pole zq follows 



oi (s - 1) + 02^0 + as (s - 1) Zo 

+ a^Zg + asZg (s — 1) + a^z^ + a-jz^ 4- Zq = 0, 



(52) 
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FIG. 4: (Color online) Fano factor F as a function of ei — £2 
calculated using the occupation-state basis {F — Fo) and the 
eigenstate basis {F — Fc) for a large interdot coupling SI = 
32 /ieV. 



with 

fli = -^rLrfl{4£272 + (71 + 72 + Tr) 

x[4f72 + 72(7i+72+rR)]}, 

02 = ^|rL(7i + 72)[4(e' + 2n^) + (71 + 72)'] 

+ {(71 +72) [4f^^ +72(71 +72)] 

+rLm^ + 3(71 + 72)2] + 4^2 + rL)}rij 

+ [4^(2 + (71 + 72)(272 + STlWI + (72 + TL)r 

03 = -rLr„[2n^ + 72(71 + 72 + Tb)1 
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7? + 72 + 4(e' + ^n^){TL + Tr) + 5TlTI 
+ 7?(372 + + Sr.^) + 7|(5rL + 7Tr) 

+72(4£2 ^ g^2 ^ lor^FR + 7r?j) 

+71 [4e2 + 8f}2 + 37I + IOTlTr + ST^ 

+i072(rL + rflj] 



(53) 



The expressions for 05, ag and 07 are not involved in 
further calculations and arc not reported here. Using also 
Eqs. P5)) and ([5T|) . the Fano factor i^o in the occupation- 



state basis is found to be 



i^o = 1 + 2 X 



On 



(54) 



Without any phonon dissipation effect, i.e., 71 = 72 = 0, 
the Fano factor becomes 



SVl'Tl [4£2 {Tr -Tl) + -iTLTl + T\ + B^^Ffl] 

[TLTl+ATLe^+m^VR + 2TL)f 

(55) 

which is identical to the previous results^^i^ obtained in 
the occupation-state basis. For a super-Poissonian noise, 
one has Fq > 1. From Eq. ([55|) . it follows that 



4e2(ri - Tr) > iVLTl + Tl + Sn'VR. (56) 

Alternatively, using the eigenstate basis, one can ob- 
tain from Eqs. ([571) ^-^d (US]) the following equation for 
the pole zq: 



where 



bi{s - 1) + b2Zo + bszo + bizo + = 0, (57) 

-rLTR[a\X2 + [^^Tr) + P^Xi], 
Xi{rL+TRf3^) + X2{TL+a^rR) 
+ rH[rL(a" + - 2a^P^) + a^P'^TR], 
~2o?(?Vj^Vr, 

Ai+A2+rL+rR. (58) 

In contrast to Eq. (j52p . only four coefficients &i (i = 1 to 
4) appear in Eq. ([57]) . From an equation for analogous 
to Eq. ([54|) for Fq, we get 



61 

bz 
64 



1 



2rir, 



{2c? (3^^ - (Ai + A2 + Fl + Ffl.) 



Y.{(?\^+c?\2^'O?0^VR\) 



(59) 



where S is given in Eq. (|41bp . Without phonon dissi- 
pation, i.e., Ai = A2 = 0, one obtains after substituting 
Eq. (USD into Eq. 



^ _ 8»^Fl X {4£2 {Vr - Fl) + %^^Vr] 
[4e2ri+4f]2(2FL-t-FR)]2 



Figure H] presents both of the calculated Fano factors 
Fo and Fe of the shot noise based on the occupation- 
state basis and eigenstate basis, respectively. At the res- 
onant tunneling point, i.e., £1 = £2, both approaches 
predict that the shot noise is sub-Poissonian. For Fo, 
in the absence of phonon-induced dissipation, a super- 
Poissonian noise can be obtained with the condition 
Eq. ([SG]) . Due to the effects of dissipation, Fo has only 
sub-Poissonian noise for the whole parameter range inves- 
tigated here. In contrast, Fe has much richer behaviors of 
super-Poissonian, sub-Poissonian, and Poissonian noise 
correlations, depending on the energy difference e\ — £2. 
Moreover, Fe, but not Fo, exhibits a double-peak struc- 
ture and an asymmetry around the dip at £1 = £2. These 



9 



1.2 



O 
^— > 
o 

o 



1.0 



0.8 



0.6 



' 1 ' 1 

0.9( 


■ 1 ■ 


\\ 


^ -1 JU.U 1 .0 


V 








X U ^ 

I.I.I. 


-200 


200 



8^ - 8^ (iLieV) 



FIG. 5: (Color online) Fano factor as a function of ei — £2 
calculated using the occupation-state basis {F — Fo) and the 
eigenstate basis {F = Fc) for a small interdot coupling fl = 
1 peV. Inset: The enlarged diagram of the Fano factor F in 
the region with |£i — £2] comparable to fi. 



features were also observed in a recent experiment (see 
Ref.i). 

The double peak in the Fano factor predicted using 
the eigenstate basis can be intuitively understood as fol- 
lows. The electrons can tunnel from the DQD to the 
right lead via two channels, namely, the ground-state 
channel and the excited-state channel. At the resonant 
tunneling point (ei = £2), the tunnel rate through the 
ground-state channel is the same as that through the 
excited-state channel. This results in a sub-Poissonian 
shot noiseji^ When ei < 62, one has Tg < Fe, and 
the electron transport through the ground-state chan- 
nel blocks that through the excited-state channel. This 
dynamical channel blockade leads to a super-Possionian 
shot noiseJ^ However, when £2 — £1 ^ ^, Fg becomes 
zero and the electron can only tunnel through the DQD 
via the excited-state channel. This single-channel tunnel- 
ing gives rise to a sub-Poissonian shot noisCfi^ as shown in 
Fig. m Similarly, when £1 > £2, one has Tg > Fg and the 
tunneling through the excited-state channel blocks that 
through the ground-state channel. The noise is super- 
Poissonian for a small energy difference £1— £2, due to the 
dynamical channel blockade. When £1 — £2 ^ fi, Fe be- 
comes zero and the electron can only tunnel through the 
DQD via the ground-state channel. This single-channel 



tunneling also gives rise to a sub-Poissonian noise. 

The asymmetry of the shot noise is caused by the re- 
laxation process induced by the electron-phonon interac- 
tion. When £1 < £2 , we have Tg < Fe and the relaxation 
process from the excited state to the ground state en- 
hances the dynamical channel blockade. However, when 
£1 > £2, one has Tg > Fe and the relaxation process from 
the excited state to the ground state suppresses the dy- 
namical channel blockade. The asymmetry of the Fano 
factor hence follows [see Fig. 

We have shown using Fig. 0] that for a large inter- 
dot coupling, e.g., O = 32 fieV, the Fano factor in the 
occupation-state basis deviates drastically from the Fano 
factor in the eigenstate basis. This verifies that the small 
interdot coupling approximation in deriving the master 
equation in the occupation-state basis is invalid. Instead, 
the master equation in the eigenstate basis should be 
used. As a further consistency check, Fig. [5] shows the 
Fano factor of the shot noise for a small interdot cou- 
pling strength (fi = 1 fieV). As expected, the calculated 
Fano factors using both basis agrees with each other ex- 
cept for the small region with |£2 — £1 1 ~ f2; in this small 
region, the results in the two cases are different because 
the condition |£2 — £i| ^ is not satisfied (see the inset 
of Fig.O. 



V. CORRECTION TERMS FOR THE 
OCCUPATION-STATE MASTER EQUATION 

In this section, we derive a controlled series expansion 
for the scattering term in the quantum master equation 
with respect to the interdot coupling strength. This pro- 
vides a concise quantitative description of the approxima- 
tion used in the occupation-state approach. In general, 
it can also allow one to derive correction terms, either 
to improve the results based on the occupation-state ap- 
proach or to estimate the resulting error. 

The master equations in both approaches arc derived 
from Eq. (|10p in the interaction picture. To study the 
difference between the two approaches, we first transform 
Eq. pIT]) back to the Schrodinger picture and get 

p<j(<) = -iTrB[i7o,Pd(tW(0)] 

00 

^Tte J dT[Hi^t,e-'""^[Hint,Pdit)pEiO)V""l, 


(61) 

where we have put t ^ t — t' and assumed t ^ 0. It can 
further be written in the more compact form^ 

Pdit) = TrE[Copdit)pE] 

00 

+TrE J dTCn,te^''^Cn,tpdit)pE, (62) 


where Cq and £int are the Liouville operators for the 
Hamiltonians Hq and TJjnt, respectively. The Liouville 
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operator /!oj for instances, is defined by CqA = — z[iJo, A] 
for any operator A. We have also used e'-'"'^A = 
^-tHoT j^^tHaT -^j^jgjj follows directly from the Baker- 
Hausdorff lemmai^ 

In the eigenstate basis, e^°'^ in Eq. (|62p is treated 
exactly. However, in the occupation-state basis, it can 
only be approximated. To illustrate this approximation 
and study the associated correction terms, we note that 
Cq ~ Ci+ Cn because of Eq. p6|) and derive the Dyson 
series 

g£or ^ ^c,r ^ j dr' e'^'^^-^'^ Cnc'^'^' + • • • , (63) 
Jo 

where Ci and Ca denote the Liouville operators for Hi 
and H^, respectively. Equation (|62p then becomes 

Pd{t) = i:vE[CoPd{t)pE\ 
Jo 

oo T 

+Tye JdrJ dT'Ci^,e^'^^~^'^Cne'''^'CintPd{t)pE 



+ ■■■ . (64) 

Taking only the first two terms, we arrive at the ap- 
proximate master equation used in the occupational-state 
basis, which is identical to Eq. (1621) with e^°'^ approxi- 
mated by e^^'^ . The third term in Eq. ([64|) is then the 
leading correction term for the master equation in the 
occupational-state basis. It consists of terms of order 
0{ililff.) or O(ilA^). Expressions for higher order correc- 
tion terms in the occupational-state approach can simi- 
larly be calculated. 

For our DQD problem, the correction terms can also be 
obtained by a direct comparison with the eigenstate basis 
result. Without the lose of generality, we assume £2 — 
El > in the following discussion. In the small interdot- 
coupling limit with f2 <C 1^2 — ct and f3 defined in 
Eq. ([55]) reduce approximately to 

a«l, /3« 77, (65) 

where r] — Cl/s <^ 1. The transformation between the 
two bases given in Eq. (|29p can be approximated by 



first order in 77, one has 

Pd{t) w -i[HuQB, Pd{t)] + —'D[a\]pd 
+— ^NT^lpd + —V[alai]p + —V[a\a2]pd 
-T LvlalPdai + a\pda2\ 

r t , t 1 71 - 72 r , 1 

-7177 [a 1 02 Prfcr^ + cr^pda^aij ^ -rWxPd + Pd<yx\ 

-72?7[4aiPdO-2 + '^zPda\a2] ■ (67) 

Indeed, when 77 = 0, this equation reduces to the ap- 
proximate master equation in the occupation-state basis 
[Eq. (|22p ]. The terms proportional to 77 are the lead- 
ing correction terms of the order ©(Jlfi^j,) or ©(fJA^) as 
expected. 



VI. CONCLUSION 

In summary, we have derived two master equations in 
both the occupation-state basis and the eigenstate ba- 
sis to describe the dynamics of the DQD. We show that 
the master equation in the occupation-state basis is only 
valid for a small interdot coupling, while the master equa- 
tion in the eigenstate basis is valid for an arbitrary in- 
terdot coupling. To demonstrate the difference between 
these two master equations, we focus on the current and 
shot-noise properties in electron tunneling through the 
DQD. When the interdot coupling is much smaller than 
the energy difference between the two dots, the current 
and shot noise in the occupation-state basis are very close 
to those in the eigenstate basis. For a large interdot cou- 
pling, however, the properties derived in the occupation- 
state basis deviate drastically from those in the eigen- 
state basis. This reveals that the master equation in the 
occupation-state basis is not accurate for the case of a 
large interdot coupling and in this case the master equa- 
tion in the eigenstate basis should be used. Also, we show 
that the shot-noise properties predicted using the eigen- 
state basis can successfully reproduce the features found 
in a recent experiment 1^ Moreover, we have discussed the 
relation between these two master equations and show ex- 
plicitly that the master equation in the occupation-state 
basis only includes low order terms with respect to the 
interdot coupling, compared with the master equation 
derived in the eigenstate basis. 



I<?>«|1>-'7|2), 

|e)«r;|l) + |2>. (66) 

Substituting Eq. into the master equation in the 

eigenstate basis [Eq. (|34p ] , and keeping terms only up to 
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Similarly, the second term in Eq. ((T4)) can be calculated 



as 



APPENDIX A: DERIVATION OF MASTER 
EQUATION IN OCCUPATION-STATE BASIS 

In this appendix, we give further details of the deriva- 
tion of the master equation in the occupation-state basis 
outlined in Sec. IIA. We first evaluate £TPd(^)- Using 
the expression for Hrit) in Eq. the first term in 

Eq. (fTil) becomes 



■^0 ^ Ik ^ 



Ik 



■ ^ nl.a^alTlT^p', {t) e*(---^)- U.Crk 



rk 



rk 



(Al) 



where t ~ t — t'. When the electron density of states in 
an electric lead is dense, each sum in Eq. (|Aip can be 
replaced by an integral. After some algebra, we obtain 



- ^ Tr„y dt'[HTmT{t')p'At) Pleads (0) 

a=Lr 

= [aialPdfi i^i) + a\aipiifi (wi) 

r r - 
— — a2alTlTrp'jr (t^2) + ala2TrTlp^Jr (wa) 

(A2) 

where a; = 2tt pi^^rk^fk rk ^'^^ electron tunneling rate 
through the left (right) barrier. Here 



(A3) 



is the Fermi-Dirac distribution with /i^ being the chem- 
ical potential of lead a and fa{'-Oi) = 1 — /q(w,;). Note 
that, in deriving Eq. (|A2p . we have used the relations 

(clfeCafe) = /a(wQfe), (Cafe^fc) = 1 - fai^ak), (A4) 



and 



/"OO 

/ ciTe±*("°'=-"''"«^(5(, 
Jo 



UJak~UJi). (A5) 







^ Tr„ I dt'[HT{t)piit) pi^^dMHTit')] 

OL—l.r 

2 



('^i) + <Apaaifi (wi) 
a2Ttp^T,.a^/r (^2) + alTrpiTla2fr {0J2) 



(A6) 



Substituting Eqs. (|J2|) and into Eq. ([13]), one ob- 
tains 



^Tp'At) - ^^P[ai]p^(t)/z(c.i) + ^V[al]p',{t)Mu,^) 



piiWfr (^2) , 
(A7) 

where V (acting on any operator A) is defined by 

V[A]p = 2ApAl - Ap - pA^ A, (A8) 

for any given operator A. 

Following similar procedures, substituting the value of 
Hcp(t) in Eq. ([2T|) into Eq. and after some algebra, 
one obtains 

>CphP^(<) = [a+] piit) + p^(0, (A9) 



with 



71 = 27r {J (e) [n (e) + 1] + J (-e) n (-e)} , 

72 - 2^ {J (-£) [71 (-£) + 1] + J (e) n (e)} , (AlO) 



where 



is the bath spectra density and 

1 

n{e) 



exp (e/kBT) — 1 



(All) 



(A12) 



is the Bose-Einstein distribution. 

With Crpiiit) and CphPdit) given by Eq. fMI and 
Eq. (|A9|) . the master equation, Eq. (fT3|) . for the reduced 
density matrix of the DQD in the interaction picture is 
found to be 



PdW = ^V[a,]pUiicJi) + '^V[a\]p'ji{u;,) 



-^V[a2TMMc.2) + ^V 
-^V[a+]pi + :^V[a^]pi 



Pdfr (^^2) 

(A13) 
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Next, we assume both a large bias voltage across the 
DQD (i.e., > uji, uj2 > ^j^b) and a very low temper- 
ature, so that fi{uJi) ~ 1, fr{^2) ~ 0. After converting 
the resulting equation into the Schrodinger picture using 
the free evolution operator e"'^"* or its approximate in 
Eq. (fT9|) . we finally have 

Pd{t) = -i[-^0-z +^Cra:, Pd{t)] 

+ llv[<j+]p,{t) + :^V[a^]p,{t), (A14) 

which is just Eq. ([22]), i.e., the master equation in the 
occupation-state basis. 

APPENDIX B: DERIVATION OF MASTER 
EQUATION IN EIGENSTATE BASIS 

This appendix gives further details on the derivation of 
the master equation in the eigcnstate basis given in Sec. 
IIB. Substituting Eq. into Eq. (Hi]), and following 

similar procedures in Sec. IIA, the dissipativc part due 
to the electric leads is evaluated to be 

^Tp'At) - ^a'V[al] pi, fiiiug) + ^f3^V[al] p^ Mco,) 

+I±a^V[a,] pi MiJ,) + ^-l3^V[a,] pi Mu;,) 

+^P^V[alTr] pi friu^,) + ^a^V[alTr] piU^e) 

+^(3^V[agTl]piU{u:,) + ^a^V[a,Tl]piULo,), 

(Bl) 

where a = cos(6'/2) and /? = sin(6'/2). In calculating 
Eq. (|B1[) . the fast oscillating terms proportional to e**^"* 
are neglected within the rotating-wave approximation. 
Similarly, from Eqs. ([55)1 and (|15p . the dissipativc part 
due to the phonon bath reads 

(0 = ^'Dia^l'] Pi{t) + ^V[<jf] piit), 

(B2) 

with the dissipation rates given by 

Ai = 2nJ{nQ)cos^ 0[n{nQ) + 1], 

A2 = 27rJ(r2o)cos^6'n(flo), (B3) 



where 

J(f7o) = 5^A2 5K-f}o), (B4) 

is the bath spectral density. 

Substituting Eqs. ((BT|1 and (|B2|l into Eq. the mas- 
ter equation for the reduced density matrix of the DQD 
in the interaction picture is 

Pdit) = ^«2l?[at] pi fiiug) + ^P^V[a\] pi Mco,) 
+ l±a'V[a,]piMLO,) + ^P'V[a,]pifi{u,) 
+ ^P^V[alTr] pi fr{oj,) + ^a^V[alTr] pi fAco,) 
+ ^(3^V[agTl] pi U^,) + ^a^V[a,Tl] pi /.(c^e), 
+ ^V[a'^^^]pi{t) + ^V[af]pi{t). (B5) 

Here we also consider the case of both a large bias volt- 
age across the DQD (i.e., p,L > ojg, uje > Pr), and a very 
low temperature, so that fi{ujg) = fi{uJe) ~ 1, fr{^g) ~ 
fr{oJe) — 0. Converting Eq. (jBSP into the Schrodinger 
picture using the free evolution operator e~'^°* with- 
out needing further approximation this time, the master 
equation of the reduced density matrix of the DQD is 
given by 

^l±.fV{a\\ pS) + ^f3^V[agTl] p,{t) 

+ ^a'V[a,Tl] pa{t) + ^P[<tL^'] pS) 

+ ^V[a^:^]pa{t), (B6) 

which is just Eq. ([34]), i.e., the master equation in the 
eigenstate basis. It should be emphasized that this mas- 
ter equation is valid for arbitrary interdot coupling, in 
contrast to the master equation in the occupation-state 
basis that is valid only for small interdot coupling. 
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